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ABSTRACT 

This  paper  presents  a  model  for  calculating  the  surface  tension  of 
pure  metals  and  their  alloys.  It  is  based  upon  the  theory  of  Eyring  et  al 
[5]  which  uses  classical  statistical  physics  to  describe  the  thermodynamic 
properties  oc  metals  in  the  liquid  state.  Calculations  show  the  surface 
tension  of  pure  aluminum  to  be  nearly  9  percent  greater  than  that  measured 
for  pure  aluminum  having  a  monolayer  offtoxide  (Al^O^) ,  is  within  about  10 
percent  of  measured  values  for  Al-xMg  and  Al-xCu,  where  x  is  the  weight 
percent  of  the  alloying  element.  -In-  the  present  calculations  0  <  x  <  8 
wt  %  for  Mg  and  0  <  x  <  30  wt  %  for  Cu  were  used.  The  values  calculated 
are  also  in  good  agreement  with  results  from  other  models.  ^The  model  was 
also  used  to  calculate  the  temperature  coefficient.  For  pure  aluminum  the 
calculated  values  fall  within  experimental  measurements,  and  exhibit  a 
slight  temperature  dependence.' 

1.  INTRODUCTION 

ACast  metal  matrix  composite  (MMC)  materials  reinforced  with  unidirect¬ 
ional  fibers  will  turn  out  to  be  a  low  cost  alternative  to  MMC  fabrication 
by  diffusion  bonding,  cold  and  hot  compaction,  and  other  methods. \  In  the 
development  of  this  kind  of  casting  technology,  one  of  the  most  important 
problems  to  be  solved  is  wetting  of  the  fibers.  From  an  experimental  point 
of  view,  wettability  is  poorly  understood.  From  the  theoretical  point  of 
view  there  are  a  number  of  different  approaches,  and  some  of  the  sophis¬ 
ticated  approaches  have  not  been  of  use  to  experimentalists.  There  are  a 
number  of  approaches  for  calculating  the  thermodynamic  and  surface  prop¬ 
erties  of  liquid  metals,  and  they  divide  into  macroscopic  and  microscopic 
models . 

♦Permanent  address:  MMCIAC,  Kaman  Tempo,  816  State  Street,  Santa  Barbara, 
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The  macroscopic  modeling  approach  is  the  most  general  and  is  based  on 
equilibrium  thermodynamics.  It  uses  the  classical  Gibbs  formulation  for  an 
interphase,  and  attempts  to  predict  interfacial  tension  or  adhesion  and  the 
surface  free  energy  of  liquid  metal  in  equilibrium  with  its  vapor,  or 
another  liquid,  or  a  solid  substrate.  The  Gibbs  model  is  an  energy  balance 
between  the  Helmholtz-free  energy  and  the  sum  of  entropy,  mechanical  work, 
surface  tension  and  the  chemical  potential  of  the  species  involved.  From 
this  balance  the  surface  free  energy  is  related  to  the  wettability  of  the 
liquid  metal  on  a  non-metallic  substrate  material,  and  wettability  is  then 
related  to  the  contact  angle  using  Young's  equation.  It  is  well  known 
experimentally  that  the  wetting  of  a  ceramic  by  a  liquid  metal  can  be 
inferred  from  changes  in  the  contact  angle  with  changes  in  the  parameters 
of  the  experiments.  Consequently,  this  king  of  modeling  is  used  to  "esti¬ 
mate"  the  wettability  of  ceramics  by  liquid  metals  under  specific  con¬ 
ditions.  What  is  most  often  the  case  is  the  Gibbs  formulation  is  used  post 
facto  to  explain  what  happened  in  an  experiment,  and,  therefore,  does  not 
really  possess  much  of  a  what  a  good  theory  should  possess:  predicting  the 
outcome  of  an  experiment! 

Various  kinds  of  micro-models  have  been  proposed  which  attempt  to 
introduce  an  atomic  or  molecular  "picture"  into  thermodynamic  arguments. 
Included  in  this  type  of  approach  is  the  notion  that  the  major  contribution 
to  adhesion  forces  arises  from  the  instantaneous  attraction  among  dipoles 
and  among  induced  dipole  interactions  between  the  liquid  and  it  substrate. 
This  method  has  been  fairly  successfully  applied  to  understand  hydrocarbons 
wetting  solid  hydrocarbons  and  similar  conditions  in  organic  chemistry. 
These  interactions  are  modeled  variously  with  the  well-known  Lennard-Jones , 
Devonshire  or  Morse  potentials,  and  with  the  London  formula  which  involves 
first  ionization  potentials  of  constituents.  Another  approach  of  the 
micro-model  type  argues  that  the  basic  parameter  behind  the  empirical 
relationship  between  surface  tension  and  heat  of  vaporization  per  unit 
atomic  surface  is  the  electron  density  at  the  boundary  of  the  Wigner-Seitz 
atomic  cell.  It  was  shown  by  Miedema  and  co-workers  [1-3]  that  the  surface 
energy  and  the  heat  of  vaporization  are  linearly  related  to  the  Wigner- 
Seitz  cell  by  empirical  scaling  rules.  Ihese  models  then  embody  a  com¬ 
bination  of  behaviors  at  the  microscopic  and  at  the  macroscopic  levels. 


Finally,  there  are  theories  of  liquid  metals  that  are  derived  ab 
initio .  These  formal  approaches  are  rigorous  mathematically  and  generally 
elegant.  Approximations  are  introduced  at  the  end,  either  in  the  form  of 
simplified  mathematics  or  over-simplified  radial  distribution  functions  and 
inter-molecular  potentials.  These  theories  are  not  really  suitable  for 
practical  approaches  to  predict  wettability  or  interfacial  energies  or  for 
guiding  experiments;  they  are  tools  to  understand  fundamental  properties  of 
liquid  metals.  The  interested  reader  is  referred  to  a  review  of  the 
literature  [4], 

The  approach  taken  in  the  work  reported  herein  is  intermediate  between 
thermodynamic  macro-models  and  detailed  theoretical  models. 

L.  THEORY 

The  theory  discussed  here  is  based  on  classical  statistical  mechanics 
first  developed  by  Eyring  and  co-workers  [5,6]  to  predict  thermodynamic 
properties  of  liquids.  Eyring' s  theory  is  founded  on  the  assumption  that 
the  metal  upon  melting  acquires  vacancies  that  are  moving  freely  through 
the  melt  and  that  there  is  short-range  order  in  the  liquid  but  no  long- 
range  order.  These  freely  moving  vacancies,  called  fluidized  vacancies, 
have  a  volume  fraction  that  is  made  to  cor  spond  to  the  volume  change  of 
the  metal  upon  melting,  or  about  3  to  4  percent  for  transition  and  non¬ 
transition  metals.  Using  statistical  mechanics  partition  functions  are 
written  which  account  for  the  gas-like  behavior  of  the  fluidized  vacancies, 
and  for  the  solid-like  behavior  of  the  liquid  metal.  The  relationship 
between  the  Helmholtz  free  energy  and  the  partition  functions  make  it 
possible  to  calculate  the  thermodynamic  properties  of  the  liquid  metal. 
Eyring  and  co-workers  have  shown  this  approach  to  be  quite  successful  in 
predicting  the  thermodynamic  properties  of  a  large  number  of  liquids, 
including  a  number  of  pure  metals.  However,  calculations  of  the  surface 
tension  of  liquid  metal  alloys  using  this  approach  has  not  been  reported. 
The  work  reported  herein  is  a  first  attempt  at  calculating  the  surface 
tension  of  pure  aluminum  and  binary  alloys  of  aluminum  such  as  Al-Cu  and 
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The  relationship  between  the  Helmholtz  free  energy  and  the  partition 
function  for  a  liquid  is 

A  =  -  kT  Inf  (I) 

where  k  is  the  Boltzmann  constant,  T  is  the  absolute  temperature,  and  f  is 
the  partition  function.  The  surface  tension  is  calculated  from 


N,v,T 


Y  =  \3ft 


where  N  is  Avogadro’s  number,  V  is  the  molar  volume  and  ft  is  the  surface 
area  occupied  by  a  monolayer  of  atoms  on  the  liquid  surface,  which  is  given 
by 

ft  =  wN  (3) 


where  id  is  the  area  occupied  by  one  atom,  and  is  the  total  number  of 

sides  available  for  atoms  on  the  liquid  surface.  To  calculate  the  Helm¬ 
holtz  free  energy,  the  partition  function  can  be  written  as  the  product 


f'  »  f  f 

B  mL 


where  f  .  is  the  partition  function  for  the  surface  atoms  and  f  is  the 
ml  d 

partition  function  for  the  bulk  liquid  atoms.  In  general,  any  partition 
function  can  be  written  as  the  product  of  partition  functions  describing 
each  process  involved  [7],  Therefore  [5] 


Nn  NrVV  N_(l-V  /V) 

f®  *  (f  )  B  S  (f  )  B  S 
B  s  g 

where  V  Is  the  molar  volume  of  the  metal  at  its  melting  point,  and 
s 


is  the  molar  volume  of  the  liquid  at  temperature  T  when  its  density  P(T)  is 
a  function  of  temperature,  and  M  is  the  atomic  weight  of  the  metal  under 


4 


consideration.  In  Equation  5,  f  is  the  partition  function  for  the  solid¬ 
like  behavior  of  the  bulk  liquid  and  f  is  the  partition  function  for  the 
gas-like  behavior.  The  partition  function  can  be  further  separated  into 
the  product  of  partition  functions,  thus 


fEinstein  frot  fvib  J(T) 


f  =  f  f  f 

g  trans  rot  vib 


showing  that  rotation,  vibration,  and  translation  motion  can  be  accounted 
for  in  this  manner.  J(T)  is  the  partition  function  for  internal  degrees  of 
freedom.  The  partition  function  for  the  oscillations  of  the  atoms  in  the 
solid  may  be  calculated  from  Einstein's  theory,  or 


E  /RT 
s 

p  _  ^ 

Einstein  ~7.  -0/T.  3 

(1-e 


whe-e  0  is  the  Einstein  temperature  [8].  For  metals,  0/T<<!  so  that 
Equation  9  reduces  to 


f_.  „  .  £exp(t  /RT) 

Einstein  s 


where  Eg  is  the  sublimation  energy  3nd  R  is  the  universal  gas  constant. 

For  pure  liquid  metals,  the  probability  of  diatomic  modules  is  ex¬ 
tremely  small,  and  those  formed  with  residual  impurities  in  the  metal  is 
also  very  small.  Consequently,  the  partition  functions  for  rotation  and 
vibration  are  unity. 

A  term  must  be  included  in  the  partition  function  given  by  Equation  7 
to  account  for  the  total  number  of  positions  available  to  an  atoms  in  the 
melt,  or  [5] 


1  +  n 


aE  V 
s  s 

(V-V  )RT 
s 
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where  n  and  a  are  parameters  that  can  be  calculated  theoretically.  It  will 
be  seen  shortly  that  in  the  further  development  of  the  theory,  these  two 
parameters  drop  out. 

The  partition  function  for  the  bulk  behavior  of  the  liquid  is  then 
written  as 


E  /RT 
s 

-e/T, 3 


(l-e  ‘) 


1  +  n 


aE  V 

e5tp  "  (V-V  )RT 


Vs/V 


2irmkT 


Nb(!-Vs/V) 


(V-V  ) J (T) 

S 


vv-v 


A  similar  reasoning  gives  the  partition  function  for  the  monolayer  or 


e;  r  /v-v 

!,  I  +  "’  — 

(l-e  \  s 


\  a’E’V 

si  s  s 

I  exP  ■  (V-V  )RT 

/  s 


N’V  /V 

S 


N'(l-Vs/V)  [rN'(V-Vs)1|l_1 


The  second  brackets  in  Equation  12  and  13  can  be  simplified  using 
Stirlings  approximation  x!  =  (x/e)  so  that 


N.  (1-V  /V) 

D  S 


N ’ (1-V  /V) 
s 


(^y 


N' (1-V  /V) 
s 


In  Equation  13,  the  primed  quantities  refer  to  the  surface  atoms  corre¬ 
sponding  to  the  un-primed  quantities.  In  the  above  derivation  vc  must  have 
conservation  of  the  number  of  atoms,  or 


N  =  N  +  N' 
B 
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where  N  is  the  number  of  atoms  in  the  bulk  and  N  the  number  of  atoms  on 
B 

the  surface  as  a  monolayer. 

Now  we  write 


Inf'  =  ln(f  f  ,)  =  In  +  In,  , 
B  mL  B  fml. 


and  using  Equations  12  through  15  gives 


V  (  f’  [1  +  n'  [(V/V  )-l]  exp[-a'E’V  /(V-V  ) RT]  ) 

■’  =  N’-f  In  -- - - - — - 5 -  +  Inf,  (1 

f  1  +  n  [(V/V  )-l]  exp [-aE  V  /(V-V  )RT ]  \ 

'S  S  S  S  S  ' 


where 


exp[E  /RT] 

-e/T  \ 3 


exp[E'  /RT] 

-0,/T,3 


3  ( 1-e  )3  =  (1-e  ~  ,x)3 

The  Helmholtz  free  energy  for  a  liquid  involving  its  surface  is 

A  =  -  kT  Inf' 


/3A 


N.v.T 
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-1  vs  3A 


V  /  \5N' 


N.v.T 


In  Equation  20,  a  random  distribution  of  vacancies  is  assumed,  that  N 
is  the  total  number  of  sites  available  for  atoms  on  the  surface,  and 


N  V 
c 

Combining  Equation  18,  19  and  20  jives  after  some  algebra  [9] 

r-.-  Crf  kT  |x.  &W  i-J 


ft)  +  ln8' 


where 


lngr  = 


1  +  n  [(V/V  )-l]  exp [-aE  V  /(V-V  ) RT ] 
s  s  s  s 

1  +  n’  [(V/V  )-l]  exp [-a ' E' V  /(V-V  )RT ] 


(20) 
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Because  n* 


=  n  and  aT  =  a,  g^=  1  and  hence  Ing^  =  0  and  Equation  22  reduces 


to 


Y  =  GJ 


-1  /Vs 


V 


kT 


E  “E'  t i  "e'/T.  f  . 

8  S  +3  In  (1~‘  1  +  lni^T) 


RT 


d-e-6/T) 


J'(T) 


(24) 


It  can  be  shown  that  ln[ J(T) /J' (T) ]  =  0.  Lu,  Jhon,  Ree  and  Eyring 

showed  that 


E' 

s 


7  E  (1  +  f) 
4  s 


(25) 


and 


0’  = 


6  (|  +  |  f)*5 


(26) 


so  that  the  second  term  in  the  bracket  of  Equation  24  reduces 


|ln^(l  +  f)  . 

Then,  we  finally  have  that 

/v  \2 

E 

y  - «  \r)  kT 

[_  4RT  0-30  +  |  1«  J(l+f)J 

In  Equation  27,  for  close  packing 


(27) 


.  -i  (vr 


(28) 


3.  NUMERICAL  RESULTS  AND  DISCUSSION 

We  will  now  perform  three  calculations:  one  for  pure  aluminum,  one 
for  Al-Cu  alloy  and  Al-Mg. 
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3.1  Pure  Aluminum 


Table  I  present  the  numerical  values  of  the  parameters  used  in  the 
calculation  of  surface  tension  of  pure  aluminum.  To  calculate  the  surface 
tension  of  pure  liquid  aluminum  the  temperature  dependent  density  values  of 
Gebhardt  et  al.  [11]  were  used.  Between  933°  and  11736K,  these  values  are 
well  represented  by 

p (T)  =  2.368  -  2.63  x  10_4(T-T  )  (29) 

m 

where  T  =  933°K  is  the  melting  temperature  of  pure  aluminum  and  T  is  the 
m 

temperature  of  interest  in  °K.  The  values  calculated  with  Equation  29 
differs  by  0.71  percent  at  the  melting  point  from  values  obtained  with  a 
similar  expression  reported  elsewhere  [12],  and  0.56  percent  at  1173°K,  the 
values  obtained  with  Equation  29  being  lower  than  those  obtained  with  a 
similar  equation  [12]  by  the  indicated  percentages.  The  molar  volume,  V  , 
for  pure  aluminum  given  in  Table  I  differs  from  that  given  by  Shimoji  [4] 
which  is  11.4  cm3/g-atom. 

Figure  1  shows  the  surface  tension  of  pure  aluminum  as  a  funct'on  of 

temperature  calculated  with  Equation  27.  These  values  are  compared  with 

other  calculated  values  [9,13.15]  and  measured  values  [9,14,16].  /Iso 

shown  is  a  range  of  values  at  660°C  of  the  surface  tension  of  pure  alum  ,mm 

measured  by  Garcia-Cordovilla  and  co-workers  [13]  using  the  method  ->f 

maximum  bubble  pressure.  They  found  that  the  surface  tension  of  pure 

_2 

aluminum  was  a  high  as  1122  mJm  [16]  for  an  unoxidized  bubble  surface, 

_2 

decreasing  to  865  mJm  with  a  controlled  increase  in  the  bubble  surface 
oxidation.  The  values  calculated  with  Equation  27,  are  in  good  agreement 
with  other  computations  using  various  modeling  approaches,  and  measured 
values  reported  in  the  literature. 

There  is  an  interesting  implication  from  the  range  of  measured  surface 
tension  values  shown  in  Figure  1  and  reported  by  Garcia-Cordovilla  et  al. 
[16].  These  suggest  that  the  values  reported  by  others  [9,14,17-26]  are 
low  because  the  surface  was  contaminated  by  oxides  during  surface  tension 
measurements.  Indeed,  measurements  of  surface  tension  of  aluminum  requires 
great  care  that  no  oxygen  is  introduced  into  the  system  except  under 
controlled  conditions. 
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The  surface  tension  value  of  1184  mJm  shown  in  Figure  1  was  cal¬ 
culated  by  Chacon  et  al  [13]  for  a  free  surface  of  pure  aluminum.  This 
value  is  thus  still  greater,  by  5.2  percent,  than  the  highest  measured 
value  [16,26]  under  conditions  of  ultra-high  vacuum,  or  with  a  clean  and 
purges  system  and  using  argon  gas  having  a  purity  of  99.9995  percent. 
Moreover,  to  attain  such  measured  high  surface  tension  values,  high  purity 

(99.999  percent)  aluminum  was  used  [26].  The  drop  in  the  surface  tension 

-2  -2 

from  1122  mJm  to  a  mean  value  of  868  mJm  [16,26]  has  been  established 
as  caused  by  a  monolayer  of  oxide,  presumably  Al^O^  [26]. 

3.2  Aluminum-Magnesium  Alloy 

For  the  case  of  binary  alloys,  Equation  27  was  modified  by  using 
rule-of-mixture  types  of  scaling.  If  the  alloy  is  designated  as  Al-xMg, 
where  x  is  the  weight  percent,  then  the  following  rules  of  mixtures  can  be 
used.  The  area  occupied  by  the  atom,  ai ,  can  now  be  assumed  to  have  the 
value 


w  =  (1-x1)oja1  +  XjUng  (30) 

where  x^  is  the  fraction  by  volume.  The  molar  volume  of  the  alloy  can  be 
written  similarly  as 

V  =  (l-x,)V  +  x, V  M  (31) 

s  1  sAl  1  sMg 

and  the  molar  volume  of  the  alloy  as  a  function  of  temperature  can  be 
written 

, .  a-*  +  *  !V_ 

1  <=A1(T>  1  °Mg<T> 

'  +  xIMMg  FTT)]  <32) 

where  is  on  the  order  of  a  few  percent  in  this  calculation.  For  the 
density  relation  we  used  the  following 


p (T)  =  2.376  -  2.8  x  10  (T-T  )  -  0.9x. 

m  i 


%  v 


where  the  last  term  Is  to  account  for  the  observation  that  the  density  of 
the  alloy  varies  linearly  with  magnesium  content  [16].  The  following 
equation  was  used  to  calculate  the  dissociation  energy  for  the  alloy  [27], 


E  =  (1-x,)  E  +  x.  E  +  2x .  ( 1-x . )  \J  E  E  .  (34) 

s  1  sAl  I  sMg  1  1  v  sAl  sMg 

These  modifications  were  made  to  Equation  27,  and  the  following  values  were 

used  in  addition  to  those  shown  in  Table  I:  E  =  32.9498  kcal/mole,  V  „ 

sMg  sMg 

=  14.821  cm3/g-atom,  M^  =  24.32  gm/mole. 

Figure  2  shows  calculated  values  of  surface  tension  for  Al-xMg,  where 
0  £  x  i  8  wt  Z  Mg.  These  results  are  compared  with  measured  values  report¬ 
ed  by  Garcia-Cordovilla  and  co-workers  for  the  unoxidized  and  the  oxidized 
state.  The  same  data  were  replotted  as  a  function  of  magnesium  content  at 
973°K,  and  are  shown  in  Figure  3.  The  results  of  the  present  calculation 
fall  between  the  results  for  unoxidized  and  oxidized  measured  values  [16], 
the  entire  range  being  about  22  percent. 

3.3  Aluminum-Copper  Alloy 

The  calculation  described  for  Al-xMg  alloy  were  repeated  for  Al-xCu 
alloy,  with  0  i  x  S  30  wt  %  copper.  The  results  of  these  calculations  are 
plotted  in  Figure  4.  These  results  are  compared  with  measurements  reported 
by  Eremenko  et  al  [28]  and  Laty  et  al  [29],  It  should  be  noted  that  the 
present  results  are  about  6  percent  higher  than  the  values  measured  at 
973°K  for  20  wtZ  copper.  Unlike  the  results  for  Al-xMg,  the  surface 
tension  of  this  alloy  increases  with  increasing  weight  percentage  of 
copper.  This  is  shown  in  Figure  5  where  the  results  of  the  present  calcu¬ 
lations  are  compared  with  the  results  of  calculations  performed  by  Poirier 
and  Speiser  [15].  In  this  calculation  the  liquidus  temperatures  of  Al-Cu 
binary  corresponding  to  the  wt%  of  copper  were  used.  These  results  give  a 
curve  which  is  parallel  to  the  curve  calculated  by  Poiriet  and  Speiser 
[15],  who  used  a  thermodynamic  model  to  obtain  their  results.  The  present 
results  give  surface  tension  values  for  Al-xCu  that  is  higher  than  Poirier 
and  Speiser’s  result  by  approximately  8  percent. 

3.4  Temperature  Coefficient. 

The  temperature  coefficient  can  be  calculated  directly  from  Equation 
27  by  dif f eriention  with  respect  to  T.  This  gives  the  following  result 


dT  =  ~  Lvcflj  ln  4  (1+f)  -  5-26  *  l0'4  (35) 

where  the  temperature  dependence  of  the  molar  volume  and  density  is  shown 
to  emphasize  the  fact  that  dy/dT  has  a  temperature  dependence.  It  is 
easily  shown  that  the  first  term  of  Equation  35  is  approximately  25  percent 
that  of  the  second  term.  Calculated  values  of  the  temperature  coefficient 
as  a  function  of  temperature  for  pure  aluminum  are  shown  in  Table  II.  From 
these  results,  we  note  that  the  temperature  coefficient  has  a  weak  depen¬ 
dence  on  temperature,  decreasing  by  5.8  percent  over  a  temperature  range  of 
300  degrees. 

It  is  worth  noting  that  the  calculated  value  of  the  temperature 
coefficient  falls  well  within  the  experimentally  determined  range,  this 
range  varying  from  a  minimum  value  of  -0.12  [24,26]  to  a  maximum  value  of 
-0.51  [18]  for  pure  aluminum,  using  various  experimental  methods.  This 
wide  range  in  measured  values  is  attributable  to  the  many  difficulties 
encountered  in  making  measurements  of  surface  tension  of  liquid  metal  even 
though  they  are  based  on  well-known  techniques  in  classical  physics. 
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.  Calculated  surface  tension  of  pure  aluminum  as  a  function  of 
temperature  using  Equation  27. 
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